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Abstract 

We present transport measurement results on superconducting Nb films with diluted triangular 
arrays (honeycomb and kagome) of holes. The patterned films have large disk-shaped interstitial 
^ . regions even when the edge-to-edge separations between nearest neighboring holes are comparable 

to the coherence length. Changes in the field interval of two consecutive minima in the field 
dependent resistance R{H) curves are observed. In the low field region, fine structures in the 
R{H) and T C {H) curves are identified in both arrays. Comparison of experimental data with 
calculation results shows that these structures observed in honeycomb and kagome hole arrays 
resemble those in wire networks with triangular and T3 symmetries, respectively. Our findings 
suggest that even in these specified periodic hole arrays with very large interstitial regions, the low 
field fine structures are determined by the connectivity of the arrays. 

PACS numbers: 74.25.F-, 74.78.Na, 74.81.Fa 
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FIG. 1: Scanning electron microscopy (SEM) image of the superconducting Nb film with a kagome 
array of holes. The center-to-center distance between nearest neighbor is 400 nm and the hole 
diameter is about 340 nm. 

I. INTRODUCTION 

Superconducting films with periodic arrays of artificial pinning sites have been extensively 
studied for a long time.- - - It is found that at the fields where the number of superconducting 
flux quantum $ = hc/2e in unit area is an integer multiple of the pinning sites, the so 
called commensurate effects such as peaks in the I C {H) and dips in the R{H) curves can be 
observed.-^ The most prevailing explanation for these effects is that the vortex lattice is 
commensurate with the underlying array and pinned efficiently at the matching fields. 

However, similar phenomena have also been observed in other systems including super- 
conducting wire networks and Josephson junction arrays.- - - In wire networks, when the 
magnetic flux through single plaquette does not equal to integer multiple of $o, supercur- 
rents must be induced along the loops to satisfy the ffuxoid quantization condition. There 
is an energy cost due to the induced supercurrents and thus results in a reduced T c . From 
this point of view, in wire networks, the resistance oscillations are caused by T c suppression 
at non-matching fields rather than enhanced pinning at the matching fields.— Interestingly, 
transitions from the pinning regime to the wire network regime can be observed in some 
superconducting film with pinning arrays. For example, with increasing hole diameter, the 
width of the strips between neighboring holes in a hole array becomes comparable to the 
coherence length at temperatures close to T c q (zero resistance transition temperature) and 
the system behaves like a wire network.-^*^ We notice that these findings are based on 
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square^^ 1 ^ and triangular— lattices. In those cases, with small edge-to-edge separation, 
the remaining geometrical structures of the patterned films, are indeed like wire networks 
with uniform width and small nodes. Thus the transition from pining array to wire network 
is straightforward. However, other series of hole arrays with large hole diameter may not 
have such direct geometrical simplification. Investigations of those hole arrays are needed 
to better understand the physics related to superconductors with micro/nano structures. 

In this work, superconducting Nb films with honeycomb and kagome arrays of holes are 
studied. The arrays can be viewed as diluted triangular arrays, for they can be constructed 
with 1/3 and 1/4 of the sites removed form the original lattice, respectively— The edge-to- 
edge separations between neighboring holes in these samples are comparable to the coherence 
length at temperature close to T c0 . Even though, unlike square and triangular arrays, the 
patterned films still have large interstitial regions. This distinct character makes the system 
different from wire networks, for in ideal wire networks, both the width of the stripes and the 
radius of the nodes are smaller than the coherence length and a uniform order parameter in 
the cross section of any stripe is expected. Large interstitial regions, on the other hand, can 
facilitate the nucleation of Abrikosov vortices^ 1 ^— , resulting in the appearance of normal 
cores. In our experiments, a series of minima are observed in the R{H) curves. While the 
oscillation period at low fields is in good agreement with the value derived from the hole 
density, the periodicity in higher fields is much larger. It is found that the transition of the 
two regions is due to the presence of interstitial vortices in the high field region. Surprisingly, 
in the low field region, wire network behaviors are observed in both samples. 

We identify the wire network behavior by the fractional matchings (fine structures) in 
the R{H) and T C (H) curves. The positions and the relative values of the matching minima 
are studied in detail. We notice that the connectivity of a wire network determines the 
characters of the fine structures, therefore we highlight the connectivity of our hole arrays 
and simplify each of them to a wire network. The comparison of the results in the holes 
arrays with those of their corresponding wire networks, including reported experimental 
data as well as calculation results based on the Alexander model,— demonstrates that the 
hole arrays studied in this work are well described by wire networks when subjected to small 
field. 
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II. EXPERIMENT 



The nano-structured superconducting films were prepared as follows. First, the super- 
conducting Nb film with a thickness of 100 nm was deposited by magnetron sputtering on 
Si substrate with SiC>2 buffer layer. Next, a micro-bridge for four terminal transport mea- 
surement was fabricated by ultraviolet photolithography followed by reactive ion etching. 
Then the desired arrays covering the whole bridge area of 60 x 60 /im 2 was patterned by 
electron-beam lithography on a polymethyl metacrylate (PMMA) resist layer. Finally, the 
pattern was transferred to the Nb film by magnetically enhanced reactive ion etching. In 
both the honeycomb and kagome samples, the value of the center-to-center distance between 
nearest neighbor (a) is 400 nm and the hole diameter (d) is about 340 nm. The scanning 
electron micrograph image of the kagome sample is shown in Fig. [TJ The smallest width of 
the stripes between the adjacent holes is about 60 nm. 

The transport measurements were carried out in a commercial Physical Properties Mea- 
surement System (PPMS) manufactured by Quantum Design. The magnetic field was ap- 
plied perpendicular to the film surface. During the measurements, the temperature stability 
was better than 2 mK. The zero field transition temperatures T c (0) are 8.713 K for the 
honeycomb sample and 8.755 K for the kagome sample, using a criterion of half the nor- 
mal state resistance at 9 K, which are 8.18 Q and 9.02 Q, respectively. The transition 
width of these two samples is about 0.15 K. The reference film without any pattern has a 
slightly higher T c of 8.87 K and transition width of 50 mK. We have measured the T C (H) 
phase boundary of the reference sample and obtained the zero-temperature coherence length 
f (0) = 9.9 nm.2 

III. RESULTS AND DISCUSSIONS 
A. Reconfiguration 

Figure [2] shows the R{H) curve measured at 8.55 K with a current of 50 fiA. A series of 
local minima and maxima are observed. Among these, seven integer matching minima can 
be clearly identified. The field interval between two consecutive integer minima versus the 
corresponding index number are shown in the inset. At fields below 300 Oe, the observed 
intervals are about 96.5 Oe, is in agreement with the value 99.6 Oe derived from the hole 
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FIG. 2: (Color online) R(H) curve of Nb film with honeycomb array of holes measured at 8.55 K 
with a current of 50 fiA. f is the filling ratio indicating the flux quanta per unit cell. The inset 
shows the field spacing of two consecutive minima as a function of the minima index. 

density. When the field is larger than the third matching field, the field spacing has a larger 
value of about 146.5 Oe. This value is very close to the matching field of a triangular lattice 
with the same lattice constant which is 149.4 Oe. One exception is that the spacing between 
the 4th and 5th minima is 100 Oe. 

Figure [3] is the R(H) curve of the kagome sample measured at 8.60 K with a current 
of 10 fiA. Several resistance minima are identified. As can be seen more clearly from the 
inset of Fig. El the field intervals have two sets of values too. The field intervals of the first 
three minima are about 108.7 Oe, again in good agreement with the first matching field 
(Hi =112.1 Oe) of the underlying kagome lattice. However, the value of the intervals jumps 
to about 148 Oe when the field is larger than the third matching field. 

The change in the field interval of two consecutive minima have been observed in rect- 
angular arrays and been explained in terms of reconfiguration of the vortex lattice from 
rectangular to square.- 1 ^ The interstitial regions in rectangular arrays are stripes which 
lead to the channeling effect of vortices.— 1 ^ While in honeycomb and kagome arrays the 
interstitial regions have disk-like shape. At the typical measurement temperatures 8.60 K, 
the diameters of the interstitial regions, approximately 520 nm, are about six times larger 
than the coherence length. Thus, Abrikosov vortices can be easily accommodated by the 
interstitial sites. 
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FIG. 3: (Color online) R{H) curve of Nb film with kagome array of holes measured at 8.60 K with 
a current of 10 [j,A. f is the filling ratio indicating the flux quanta per unit cell. The inset shows 
the field spacing of two consecutive minima as a function of the minima index. 

Form the field spacing values between two consecutive integer matchings, the overall 
feature of the R{H) curves can be explained as follows. When the field is small, supercurrents 
are generated around the holes in order to satisfy the fluxoid quantization condition. At 
the integer matching field, it is the arrangement of the fluxoids that is commensurate with 
the hole lattice. When the applied magnetic field is larger than the third matching field, 
Abrikosov vortices are generated in the interstitial region to minimize the free energy— These 
interstitial vortices are effectively pinned by the repulsive interactions from the fluxoids in 
the holes, which create a caging potential.— 1 ^ The overall flux lattice, which is composed of 
interstitial vortices and the fluxoids in the holes, would have a triangular symmetry. That 
is why the field intervals of the minima observed in high field region in both lattices are in 
good agreement with the value derived form a vortex lattice with triangular arrangement. 
The saturation number n s , namely, the maximum number of fluxoids trapped by one hole 
without the entry of interstitial vortices are three in both samples. This value exceeds the 
limit n s = dj (4£(T)) ~ 1 given by Mkrtchyan and Schmidt,— where d is the diameter of the 
defect and £ (T) is the temperature dependent coherence length. This fact indicates that the 
proximity of the holes is prominent.— 

However, the exact vortex configurations and their evolution are strongly dependent on 
the density of flux quanta and the confinement of the mesoscopic structure. Both im- 
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age experiment^ and simulatior*2ii2&2£ ]2 ave revealed the complex nature and the transition 
process between various configurations of the composite vortex lattice. This may also re- 
sponsible for the observation of a field interval of 100 Oe between the 4th and 5th minima 
in Fig. H 

B. Fractional matching and wire network behavior 

In the low field region, fine and repeatable sub-minima are observed as can be seen from 
Fig. [2] and Fig. [3J Because no interstitial vortex is involved in this region and the proximity 
of the holes, we can related the systems to wire networks despite the large disk-shaped 
interstitial sites. 

We notice that a wire network is obtained by assigning nodes to the center of the in- 
terstitial regions in the original hole array and connecting them. In Fig. HJ the centers of 
the interstitial regions are regarded as nodes. Then, the superconducting stripes between 
neighboring holes are viewed as wires (dotted lines) connecting the nodes. Following this 
procedure, the connectivity of the patterned films are highlighted and the honeycomb and 
kagome hole arrays are transformed to triangular and T3 wire networks, respectively. Al- 
though C. C. Abilio et al. have found that superconducting Al film with a square array of 
holes can be described as a square wire network,— the simplification we made here is more 
radical for the diameter of the disk-shaped interstitial regions is even larger than that of 
the holes. Assuming that the distance between the centers of the neighboring holes is a in 
the hole array, then in the dual lattices the side length I of the elementary triangles will be 
\/3a and the side length of the rhombus tile in the T 3 geometry will be 2\/3a/3 as can been 
seen from Fig. Ufa) and Fig. Hfb), respectively. If the matching effects of the hole arrays 
originate from T c suppression at non-matching fields, the experimental results should show 
similar features with what have been observed in their corresponding wire networks. 

Figure G^a) shows the fine structures of the R{H) curve of the honeycomb sample mea- 
sured at 8.60 K with a current of 10 /xA. The x axis is given as filling ratio, / = $/$o ; 
where $ = ifa 2 3\/3/4 is the magnetic flux per elementary triangle of the corresponding 
wire network. In the field range < / < 1, resistance minima are observed at / = 1/4, 1/3, 
1/2, 2/3 and 3/4. The magnitude of the oscillation is about eight percent of R^. At the 
field / =1/2 where the commensurate effect is most pronounced, the resistance drops to a 
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FIG. 4: Illustration of the transformation from hole arrays (dark circles) to wire networks (dotted 
lines), (a) Honeycomb hole array to triangular wire network, (b) Kagome hole array to T3 wire 
network, a: the distance between the centers of the neighboring holes. I: the side length of the 
corresponding wire networks. 

value comparable to that at zero field. At / = 1/4 and 3/4, the resistances are about three 
percent of Rjy. 

Figure [5^b) shows the T C (H) curves in reduced units, AT c /T c (0) vs /, using different 
resistance criteria r = R/R^. The curves were obtained with a current of 10 fiA. The 
parabolic background-^ which reflects the contribution of the finite width of the strands to 
the critical field has been subtracted. For r =0.05 and 0.1, dips are observed at filling ratios 
/ = 1/4, 1/3, 1/2, 2/3 and 3/4. With the criterion of r = 0.5, only dips at / = 1/4, 1/2 
and 3/4 are visible and those at / = 1/3 and 2/3 are missing. For r = 0.8, T c changes 
smoothly with field and no features can be identified. Thus, comparing Fig. [3^a) and (b), 
one can see that minima in the R{H) curve is located at the same fields where dips are 
observed in the T C (H) curve. The fine structures are more pronounced in R{H) than in 
AT c /T c (0) vs /. This is because the resistance is measured at temperatures slightly above 
T c0 in the superconducting transition region where resistance drops sharply with decreasing 
temperature, a small change in T c will result in an enhanced effect in resistance, and thus a 
more pronounced commensurate effect. 

To our knowledge, fractional matching effects in honeycomb arrays have only been ob- 
served at / = 1/2 in previous studies.— 1 ^ At the typical measuring temperature T = 8.6K 
(0.987 T c ) in this work, f (T) = £(0)/y/l - T/T c = 86.9 nm is slightly larger than the width 
of the narrow stripes. This fact makes the samples investigated here be qualitatively differ- 
ent from those in previous works on hole arrays and be better described as wire networks. 
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FIG. 5: (Color online) (a)Low field R{H) curve for the sample with a honeycomb hole array 
measured at 8.60 K and / =10 /j,A. The field is normalized by the first matching field and the 
resistance is normalized by Rn at 9 K. (b)The field dependent transition temperature T C (H) 
curves of the honeycomb sample with different criteria r = R/R^. The parabolic background 
has been subtracted. The top one is the theoretical curve for triangular wire network in unit of 

i 2 /m 2 - 

Rich structures in the T C (H) curves of wire networks may serve as fingerprints to differenti- 
ate one array from another, since for a given geometry the fine structures are only expected 
at a particular series of filling ratios.-^^ From the linearized Ginzburg-Landau equations, 
Alexander had derived the relation of the order parameters at the nodes subjected to non- 
integer flux.— Based on those equations, the task of finding the field dependent transition 
temperature is reduced to eigenvalue problems. The mathematical treatment is very similar 
to that of the tight bounding electrons in 2D arrays subjected to an external field which leads 
to the famous Hofstadter butterfly energy spectrum.— The top curve plotted in Fig. \5HJd) is 
the theoretical values of T C (H) of the triangular wire network, which is the dual structure 
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of honeycomb hole array. The curve is calculated by 



T c (o)-T c (f) m 2 
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where / is the side length of the wires, e t is the eigenvalue^ of the following equation obtained 
by using Landau gauge and periodic boundary conditions, 



it 



= 2cos[7r(2n - 1)/ - k y /2}e 
+ 2 cos [47m/ — k y ] ip n 
+ 2cos[7r(2n + l)/-/c ?/ /2]e' 




'n-l 




n+1, 



(2) 



where k y = 2%^-, k = 1,...,N implies the periodic condition, n denotes the node index 
and ip n is the order parameter at node n. The fine structures of the T C (H) of triangular wire 
network have also been studied by analytical approach based on multiple-loop Aharonov- 
Bohm Feynman path integrals.— The main features in T C (H) such as the position and relative 
strength of the most prominent dips are the same as those in the Alexander's treatment. 
The most pronounced dips occur at / = 1/4, 1/3, 1/2 etc., in good agreement with the 
experimental data. At the fields where dips are observed, the order parameter at the nodes 
interferences constructively and form different locked-in states corresponding to local max- 
ima in T c .— In the fabrication process, contamination and damage of the samples are most 
significant to the narrow strips, resulting in a lower T c and a widening in the transition 
width compared to the reference film. The phase coherence between adjacent nodes is fully 
established when all the stripes are in superconducting state. Therefore, the interference 
effects is strong near the zero resistance state r = 0.05 and significantly reduced for r > 0.1. 

Figure |6](a) is the low field R{H) curve of the kagome sample measured at 8.67 K. With 
kagome geometry, / corresponds to the magnetic flux through a rhombic tile with side length 
2\/3a/3: / = $/$ — 2y / 3a 2 i7/(3 < ! ) o)- Again we focus on the results for / in the interval 
between and 1. Dips are identified at /=l/6, 2/9, 7/9 and 5/6. These fine structures are 
observed for the first time in kagome hole array, perhaps because of the specified geometrical 
parameters. Anomalies are also visible at / = 1/3 and 2/3, but they are more like kinks 
than dips. Most strikingly, in contrast to the dip observed in the honeycomb sample, a 
remarkable peak is observed at / = 1/2. The resistance is far above the background, 
reaching 16.9 percent of the normal state resistance. 
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FIG. 6: (Color online) (a)Low field R(H) curve for the sample with a kagome array measured at 
8.67 K and 7=10 /xA. (b) The field dependent transition temperature T C (H) of the kagome sample 



with different criteria r = R/R^. The parabolic 



the theoretical curve for T3 wire network (Ref. 



background has been subtracted. The top one is 



3a) in unit of Z 2 /£(0) 2 . 



Figure [6]^b) shows the T C (H) curves in reduced units, AT c /T c (0) vs /, of the kagome 
sample determined by using different criteria of r. Again, all the measurement were carried 
out with a current of 10 //A. For small value of r (0.05 and 0.1), dips present at 1/6, 2/9, 
1/3, 2/3 and 5/6, although the one at 2/3 is relatively weak (see the labels in Fig. E](b)). 
For larger values of r (0.5 and 0.8), the dips at 1/6, 2/9 and 1/3 gradually disappear but 
the peak at / = 1/2 remains pronounced. 

Most of the fractional matchings observed here are absent in previous works performed 
on kagome lattice.— In a recent work on hole arrays with small hole size which is 
in the pinning regime limit, clear fractional matchings were only observed at / = 1/3 and 
2/3.— Contrasting to their findings, in our work with edge-to-edge separation comparable 
to the coherence length, dips at / = 1/6 and 5/6 are most pronounced and a distinct peak 
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anomaly at / = 1/2 is observed. The results agree well with what have been observed in T3 
wire networks,—^ the dual structure of the kagome hole array as seen in Fig. H](b). The top 
curve of Fig. El^b) shows the theoretical curve of T 3 networks.— This is done by the following 
equation which relates the eigenvalues of T 3 (e(/)) at filling ratio / to the eigenvalues for 
triangular lattice at 3//2.— 

£ 2 (/)-6 = 2cos(7r/)e t (3//2) (3) 

The comparison with the theoretical curve shows that only the dips at f=7/9 and other 
weaker ones are absent in the experimental curves. For /=l/6, 2/9 and 1/3, constructive 
interference occurs among the superconducting order parameters of the interstitial sites 
and the fluxoids establish long range locked-in commensurate order. While for f=l/2, 
superconductivity is localized in single tiles and long range coherence between network sites 
cannot be established.— This kind of fully destructive quantum interference has also been 
observed in kagome wire networks.— 1 ^ The oscillation of R{H) curve with a peak at half a 
flux per tile is similar to the single- loop Little-Parks effect. In that case, the supercurrent 
density reaches the largest value at half flux in order to satisfy the fluxoid quantization 
condition.— 

The results obtained from these two lattices confirms the validity of the transform from a 
hole array to a wire network. However, in contrast to the small nodes in ideal wire networks 
or square arrays of holes, the geometries studied here possess very large superconducting 
disk-shaped nodes which can trap Abrikosov vortices easily. The formation of Abrikosov 
vortices implies spatial variation and presence of zero points of the order parameters in the 
interstitial region. Then the simplification of the large interstitial region to a node in a 
wire network is no longer appropriate. Therefore, transformation from hole arrays to wire 
networks is only valid for small field values. 

IV. CONCLUSIONS 

In conclusion, we have studied the commensurate effects in superconducting films with 
honeycomb and kagome arrays of holes with small edge-to-edge separation. We found that 
the magnetoresistance curves have two regions with different oscillation periods. At small 
fields, there are no Abrikosov vortices presents in the interstitial regions. The large disk- 
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shaped interstitial region can be simplified as a single node and a one-to-one correspondence 
between hole arrays and wire networks is established. Comparison of experimental data with 
calculation shows that the simplification works well. Our results suggest that at low fields, 
the behavior of these specified periodic hole arrays are determined by the connectivity of 
the systems. 
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